hep-ph/0703092 
HRI-P-07-03-001 
OUTP 0703P 



Turbulent Supernova Shock Waves and the Sterile Neutrino Signature 

in Megaton Water Detectors 

Sandhya Choubey 1 *, N. P. Harries 2t , G.G. Ross 21 

1 Harish- Chandra Research Institute, 
Chhatnag Road, Jhunsi, Allahabad 211 019, INDIA 

2 The Rudolf Peierls Centre for Theoretical Physics, 
University of Oxford, 1 Keble Road, Oxford, 0X1 3NP, UK 



Abstract 

The signatures of sterile neutrinos in the supernova neutrino signal in megaton water Cerenkov detectors 
are studied. Time dependent modulation of the neutrino signal emerging from the sharp changes in the 
oscillation probability due to shock waves is shown to be a smoking gun for the existence of sterile neutrinos. 
These modulations and indeed the entire neutrino oscillation signal is found to be different for the case 
with just three active neutrinos and the cases where there are additional sterile species mixed with the 
active neutrinos. The effect of turbulence is taken into account and it is found that the effect of the shock 
waves, while modifed, remain significant and measurable. Supernova neutrino signals in water detectors can 
therefore give unambiguous proof for the existence of sterile neutrinos, the sensitivity extending beyond that 
for terrestial neutrino experiments. In addition the time dependent modulations in the signal due to shock 
waves can be used to trace the evolution of the shock wave inside the supernova. 
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1 Introduction 



The detection of neutrinos from SN1987A in the Large Megallenic Cloud remains a major landmark in neutrino 
physics and astrophysics. The 11 events detected in Kamiokande 1 and 8 events in 1MB [2] stimulated a 
plethora of research papers exploring both type-II supernova dynamics and neutrino properties [3] . A supernova 
explosion within our own galaxy will generate tens of thousands of events in the currently running and proposed 
neutrino detectors and hence is expected to significantly improve our understanding of the type-II supernova 
explosion mechanism on one hand and neutrino physics on the other. 

Our knowledge on the pattern of neutrino mass and mixing has seen tremendous improvement from the 
results of a series of outstanding experiments with solar [4], reactor [5j [6], atmospheric [7] and accelerator 
8, 9 J neutrinos. The existence of neutrino flavor mixing and oscillations have been established beyond doubt. 
The 3cr range for the mixing parameters governing solar neutrino oscillations is Am\ x — (7.2 — 9.2) x 10~ 5 
eV 2 and sin 2 #12 = 0.25 — 0.39 [10j . while the atmospheric neutrino oscillation parameters are restricted to be 
within the range Am\ x = 2.0 - 3.2 x 10~ 3 and 0.34 < sm 2 6> 23 < 0.68 with sin 2 26> 23 < 0.9 [ID]. However, 
the last mixing angle #13 remains unmeasured, although we do know that it must be small. The current 3er 
upper bound on the allowed values for this parameter is sin 2 613 < 0.044 [10]. If #13 is close to the upper 
bound it opens up the possibility of observing CP violation in the lepton sector - the CP phase Sep being the 
second missing link in our measurement of the neutrino mass matrix. Also still unknown is sgn(Am§ 1 ), which 
determines the ordering of the neutrino mass spectrum, aka the neutrino mass hierarchy. The case Am 2 ^ > 
corresponds to the normal mass hierarchy (N3) while Am^ < corresponds to the inverted mass hierarchy 
(13) H A number of suggestions have been put forward to determine these three unmeasured parameters. These 
include the measurement of antineutrinos from reactors with detector set-up aiming to achieve sub-percent level 
in systematic uncertainties [IT], using intense conventional beams produced by decay of accelerated pions 
[T2I [TB] , using v e beams produced by decay of accelerated radioactive ions stored in rings ( "Beta-Beams" ) [H] 
and using v e and beams produced by decay of accelerated muons stored in rings ( "Neutrino Factory" ) [T5] . 
Most of these proposed experiments are expected to be very expensive as well as technologically challenging. 

In principle, information on two of the three parameters mentioned above can be obtained by observing the 
neutrinos released from a galactic supernova. Detailed studies on the potential of using supernova neutrinos 
to unravel the mass hierarchy and #13 have been performed in the context of water Cerenkov detectors in 
[TP] and (large) scintillator detectors in [17j . It was shown that the neutrino telescope, IceCube, even though 
designed to observe ultra high energy neutrinos, could be used very effectively to detect low energy supernova 
neutrinos [IHlQjj]. Most studies on determining of sgn(Am| 1 ) and #13 with supernova neutrinos depend on 
the hierarchy between the average energies of u ei v e and v x , where v x stands for v^, D^, is T , or v T . However 
the most recent supernova models which take into account the full range of the significant neutrino transport 
processes, predict very little difference between the average energies of v e and v x [2D]. At the same time they 
also seem to be inconsistent with equi-partition of energy between the different neutrino species, an assumption 
of all earlier papers on supernova neutrinos [21] , In a nutshell, model uncertainties in the supernova parameters 
could washout the oscillation effects and render this method of hierarchy and 6*13 determination useless. 

Less model dependent signatures of sgn(Am| 1 ) and #13 can be seen through the "shock effects" in the 
supernova neutrino signal in large detectors [TDJ [22] [23] [23] [25] [2D] . The steep density profile at the shock 
front results in a drastic change in the oscillation probability and hence can yield information on the neutrino 
properties {sgn^Am^) and #13) despite the uncertainties associated with the supernova dynamics and neutrino 
transport. However it is also expected that the shock, as it moves outward, will leave behind a turbulence in the 
density profile of the supernova matter. Early studies indicated that this could severely obscure the oscillation 
signals due to the shock wave [27] . Here we reexamine this in detail following a recent analysis of the the effect 
of turbulence in [28j . In this the nature of the shock plays an important role. In addition to the initial forward 
shock it is expected from supernova simulations that a reverse shock is formed [25] . We find that the effects of 
the reverse shock are wiped out by turbulence but that the effects of the forward shock, while changed, are still 
significant and leave a clear signal of the resonant oscillation. 

Another important question in neutrino physics is how many (if any) sterile neutrino species there are. The 

1 We adopt the convention in which m?. = m? — m?. 
Note that the same notation is used even when the neutrino mass spectrum is quasi-degenerate since the NH and IH terms 
here refer only to the sgn(Am| 1 ). 
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first, and so far the only, experimental evidence which requires the presence of sterile neutrinos comes from the 
LSND experiment [25] , While the so-called 2+2 and 3+1 scenarios involving only 1 sterile neutrino [3D] are 
now comprehensively disfavored by the global neutrino data, the 3+2 scheme with 2 sterile neutrinos mildly 
mixed with the 3 active neutrinos has a more acceptable fit to all data including LSND [31] although the value 
of the associated LSND mixing angle is still problematic [32] . The LSND result will be checked by the on-going 
MiniBOONE experiment [33] which is expected to give the final verdict at least on the oscillation interpretation 
of the LSND data. However, even though MiniBOONE should refute the LSND result, it will still leave ample 
room for the existence of extra sterile neutrinos with mixing angles below its sensitivity limit. We show that 
measurement of supernovae oscillation effects provide a much more sensitive probe for sterile neutrinos and 
analyze in detail the signal that results assuming the 3+2 LSND scenario. 

In [19] we analyzed the signature in the IceCube detector of neutrino oscillation occuring within a galactic 
supernova, taking shock effects into account and assuming that the 3+2 mixing scheme was true. We considered 
all possible neutrino mass spectrum with 3 active and 2 sterile neutrinos and showed that IceCube could 
distinguish between these different mass spectra. We also showed how well IceCube could distinguish the case 
where there were only 3 active neutrinos from the ones where there are 2 additional sterile neutrinos. In 
this paper we repeat this analysis for the case in megaton water Cerenkov detectors, including the effect of 
turbulence. Water Cerenkov detectors have a good energy resolution and allow for the reconstruction of the 
energy of the incoming neutrinos. Thus, for such detectors, we will have information on both the time as well 
as the energy of the arriving supernova neutrinos. We calculate the number of supernova v e events in a generic 
megaton water Cerenkov detector. We calculate this expected signal first for the case of three active neutrinos 
and then for the case with two extra sterile species. We show that the presence of sterile neutrinos changes the 
time evolution of the average energy as well as the flux of the supernova neutrinos. Therefore, this can be used 
to indicate the presence of sterile neutrino and also to distinguish between the different possible mass spectrum 
scenarios. The effect of the shock is shown to be significant on the resultant signal. We discuss the issue of 
turbulence and take into account the turbulence due to the passage of shock. 

The paper is organized as follows. We begin in Section 2 with a brief review of neutrino oscillations inside 
supernova. We discuss the effect of the shock wave on the oscillation probability, with and without the turbulence 
caused by the passage of shock. The impact on the neutrino oscillation probability by the turbulence in the 
matter density behind the shock is clearly outlined. In Section 3 we discuss the neutrino spectrum produced 
inside the supernova and their mode of detection in water Cerenkov experiments. In section 4 we present our 
results and discuss the implications in Section 5. We end with our conclusions in Section 6. 



2 Neutrino Oscillations Inside the Supernova 

When neutrinos propagate in matter, they pick up an extra potential energy induced by their charged current 
and neutral current interactions with the ambient matter [34|, I35|, 136] . Since normal matter only contains 
electrons, only v e and v e undergo both charged as well as neutral current interactions, while the other 4 active 
neutrino species have only neutral current interactions. The evolution equation of (anti)neutrinos inside the 
supernova can be written in the flavor basis as 

n = ^(UM 2 U^ +A) , (1) 

where U is a unitary matrix and is defined by \i>i) — Ui a \u a ), i>i and v a being the mass and flavor eigenstates 
respectively. For the case where we consider two extra sterile neutrinos, i = 1 — 5 and a—e, /i, r, s± or s 2 , where 
si and S2 are the two sterile neutrinos. The matrices M 2 and A are given respectively by 

M 2 = Diag{m\,m\,m\,m\,m\) , (2) 

A = Diag(A l ,0,0,A 2 ,A 2 ), (3) 
Ax = A cc = ±V2G FP N A Y e x 2E , (4) 
A 2 = -A NC = ±V2G F pN A (l — Y e ) x E . (5) 
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The quantities Acc and Anc a re the matter induced charged current and neutral current potentials respectively 
and depend on the Fermi constant Gf, matter density p, Avagadro number Na, electron fraction Y e and energy 
of the neutrino E. The "+" ("— ") sign in Eqs. (|4]) and ([5]) corresponds to neutrinos (antineutrinos). 

The corresponding expressions for the three neutrino framework follows simply by dropping the extra 2 
states due to the sterile components. Note that we have recast the matter induced mass matrix such that the 
neutral current part AjyCj which appears for all the three active flavors, is filtered out from the first three 
diagonal terms and hence it appears as a negative matter potential for the sterile states which do not have any 
weak interactions. Therefore, for the three active neutrino set-up the mass matrix is 

A4 2 = Diag{m\ 1 m 2 , m^) , (6) 

A = Diag(A 1 ,0,0), (7) 

The mixing matrix U is given in terms of mixing angles and CP-violating phases. If CP conservation is assumed 
the mixing matrix takes the form 

n n — 1 

u = n n r ab{9ab) , (8) 

B>AA=1 

where n = 3 for only three active neutrinos and n = 5 when two additional sterile neutrinos are present and 
R AB is an n x n rotation matrix about the AB plane. 

Assuming that all phases get averaged out, the survival probability P ee of v e after they have propagated 
through the supernova matter can be written a^l 

Pee = Y,\ U ^\ 2 \ U ™\ 2p ^ ( 9 ) 

where, Pij — |(£'j|£'£"')| and U e i and C/™ are the elements of the mixing matrix in vacuum and inside matter 
respectively and K^/l^™)! 2 is the effective "level crossing" probability that an antineutrino state created as 
inside the supernova core emerges as the state \i>i) in vacuum. Largest flavor conversions occur at the resonance 
densities. For two flavor oscillations, the resonance condition for antineutrinos involving only either the v e or 
the sterile states is given by [35] 

\A k \ = (-l) fc Am^cos2% , (10) 

where k could be either 1 or 2, Am^ = to 2 — m 2 is the mass squared difference and Oij is the mixing angle 
between the two states in vacuum. Note in Eq. (|10|). we included the sign factor (— l) fe in order to take into 
account the fact that Acc < and Anc > for the antineutrinos. Since Acc < (Anc > 0), the condition 
of resonance is satisfied only when the relevant Am^ < (Am^ > 0). For the resonance between the v e and 
sterile states one must satisfy the condition 

A C c ~ A NC = Am^ cos 2% , (11) 

which is satisfied for antineutrinos when Am^ < since Acc > Anc 

The probability of level crossing from one mass eigenstate to another is predominantly non-zero only at 
the resonance. In the approximation that the individual two flavor resonances are far apart, the effective 
level crossing probability can be written in terms of the individual two flavor "flip probability" Pij. The "flip 
probability" between the two mass eigenstates used in this paper is given by |38j 

_ exp(-7sin 2 %) - exp(~7) 
1 - exp(- 7 ) 



2 

7 = 7T- 



Am 



E 



dlnA 



dr 



(13) 



3 Since we are interested in the detection of supernova neutrino in water Cerenkov detectors where the largest number of events 
are expected through the capture of P e by protons, we will mostly refer to the survival probability of u e (Pee) in this section. 
However, similar expressions are valid even for the v e survival probability. 
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where A is Aqc-, j4jvc ot Aqc — A^c depending on the states involved in the resonance, r mva is the position 
of the maximum violation of adiabaticity (mva) [37] and is defined as 

A(r mva ) = Am% . (14) 

One can see from the above equations that Pj,- and hence the transition probability depends crucially both on 
the mixing angle as well as on the density gradient. If for given Amjj and 8ij, the density gradient is small 
enough so that 7 ^> 1 and 7sin 2 # ^> 1, P^ ~ and the transition is called adiabatic. On the other hand, 
if the density gradient is very large or Oij is very small, then Py ~ cos 2 6y . In this case we have an extreme 
non-adiabatic transition. For all intermediate values of the density gradient and Py ranges between [0-1]. 

2.1 Neutrino Oscillations with Static Density Profile 

Three active neutrinos only: As discussed above, in the three neutrino framework we have two possibilities 
for the neutrino mass spectrum, the normal (Am^ > 0) and inverted (Am^ < 0) hierarchy. We will call them 
N3 and 13 cases respectively. Also, since here we have 2 independent AmJ,, we can have 2 resonances and hence 
the final level crossing probability will be given in terms of the 2 individual flip probabilities. However, it is 
now known at more than 6a C.L. from the solar neutrino data that Am^ > 0, and hence the Am^ driven 
resonance condition given by Eq. (flQ|) is satisfied only for the neutrino channel. But the sign of Amjj is still 
unknown. Hence the Ato 2 ^ driven resonance can appear in either the neutrino or antineutrino sector depending 
on whether the mass hierarchy turns out to be normal or inverted respectively. Therefore, for Ato 2 ^ < 0, the 
u e survival probability (neglecting earth matter effects) is given by 

Pee =Pl 3 |C4l| 2 + (l- Pl3)\U e3 \ 2 (15) 

where Pi 3 is the flip probability between the 1 and 3 states and can be calculated using Eq. (fT2"]) with Am 2 ^ = 
Am|i and 6ij = #13. For Ato 2 ^ > the resonance does not occur in the antineutrino channel and therefore 

Pee = \Uel\ 2 (16) 

Three active plus two sterile neutrinos: For the three active and two sterile neutrino framework we have 
4 independent mass squared differences and hence can have as many as six possibilities for the mass spectrum 
which we call [19] : 
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with Am|i > always. Since we have 4 different independent Am 2 ; and both active and sterile components, we 
can have many more resonances in this case. In fact, we expect to have resonances between the 1-2, 1-3, 1-4 and 
1-5 states due to the charged current term and between 2-4, 2-5, 3-4 and 3-5 states due to the neutral current 
term. Again since Am|i > 0, the 1-2 resonance happens only for neutrinos, while the other resonances can 
happen in either the neutrino or antineutrino channel depending on whether the corresponding mass squared 
difference is positive or negative. Another difference between the only active and active plus sterile cases is 
that while for only active neutrino cases it was enough to compute just the survival probability P ee , with sterile 
neutrinos in the fray, the total neutrino flux produced inside the supernova does not remain constant as the 
active flavors also disappear into sterile species. Thus we will need to calculate the probabilities P ee as well as 
P xe . The oscillation probability in general is given by 

P a0 = J2 P *i P %' ( 23 ) 
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Table 1: The probabilities P™ and P™ for three active neutrinos plus two sterile neutrinos, where P™ is given 
by Eq. (|25|) and Py is the flip probability at the resonance between the Vi and Vj mass eigenstates. 
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where and P^ are respectively the v a —> u™ transition probability in the supernova and v~ — > vp transition 
probability in the earth. In the absence of earth matter effects 

P?e=\U ei \ 2 (24) 

and 

3 

where Ptj is the flip probability between the i and j mass eigenstates, given by Eq. (fT2|). We explicitly give the 
probabilities PJ™ and P£j for all the possible mass spectra in Table [TJ For further details and the level crossing 
diagrams for each of these mass spectra, we refer the reader to [19] . 

2.2 Effect of the Shock Wave on P ee 

It is believed that when the core of a collapsing supernova reaches nuclear density the collapse rebounds forming 
a strong outward shock. The shock stalls and is eventually regenerated via neutrino heating. Numerical 
simulations show that as well as a forward shock a reverse shock forms and the region behind the shock wave is 
highly turbulent |25j . The detailed density profiles from numerical simulations are not available to us, therefore 
we consider the simplified profile for the shock wave used in [23j [26] . A schematic diagram showing the density 
change at the shock front for forward shock and for forward and reverse shock is shown by the dashed black 
lines in Figs. [Ha) and[^a). 

As can be seen from these figures, the effect of the supernova shock wave is to generate very sharp changes 
in the density gradient, which is expected to change the flip probability. Since, as mentioned above, the flip 
probability is calculated at the position of resonance, the effect of the shock will be to modify the flip probability 
and hence P ee when the shock crosses the resonance density for a certain Am^. Also, we see from the figures 
that the shock generates fluctuations in the density gradient. This results in the same Am^ producing multiple 
resonances which are relatively close together. If we assume that the phase effects can be neglected even in this 
case (39], then the individual resonances can be considered as independent two generation resonances and the 
net flip probability Pjj can be expressed in terms of the multiple flip probabilities Pi as (40] [23] 

(1-P H Ph \ tt (I - Pi Pi 
\ P H I-Ph 11 V Pi I- Pi 

where n is the number of resonances occurring for the same Am^ due to the shock effect. In particular, if the 
mixing angles are sufficiently large then in the region of the supernova where the density can be approximated as 
almost static, the flip probability between the mass eigenstates is approximately zero (the resonance is adiabatic). 
As the shock wave passes through the resonant density, the flip probability increases to approximately one (the 
resonance becomes strongly non- adiabatic). As the shock passes over, then in the approximation that there is 
no turbulence (which will be considered in the next subsection), the flip probability goes back to zero. The 
resultant flip probability has a time dependence shown in Figs. [IJb) andOJb). 

Since the resonance densities are different for the different Am]j involved, the shock wave will cross them 
at different times. Therefore, one can follow the evolution of the shock wave by following the time profile of the 
supernova signal [19j . Also, since the resonance density is determined by the energy of the neutrino, different 
energy neutrinos have their resonance position as different points and thus get affected by the shock at different 
times. Thus the shock imprints its signature on both the energy spectrum as well as the time profile of the 
neutrinos arriving on earth. The effect of the shock wave is discussed in [22] [23] for the forward shock and 
in [25| [26] when there is a reverse shock as well. All these papers used large/megaton water detectors for the 
supernova neutrino signal on earth. The possibility of using shock waves to discern the existence of eV mass 
sterile neutrinos in the supernova signal in IceCube was discussed in [T9"] . 

2.3 Turbulence 

Numerical simulations have shown that a highly turbulent region forms behind the shock wave [41] . The 
fluctuations in these turbulence cause flips between the mass eigenstates [371 12H]- The length scale of the 
fluctuations for which the flips are dominant are at the length scale 
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Figure 1: (a) The density profile of the forward shock as a function of the distance from the core of the supernova, 
(b) The resulting flip probability as a function of time. The dashed black line is for the case of no turbulence 
and the solid red line includes the effect of turbulence. 
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Figure 2: (a) The density profile of the forward and reverse shock as a function of the distance from the core of 
the supernova, (b) The resulting flip probability as a function of time. The dashed black line is for the case of 
no turbulence and the solid red line includes the effect of turbulence. 



Am 2 sin 29 

For the atmospheric resonance Lji uct = 47km 15 ^ IeV s i n ° 2 e 13 > which is smaller than the resolution of the 
numerical simulations and therefore a physical model of the fluctuations is required. Analytical solutions for a 
density fluctuating medium have only been found for the un-physical delta-correlated noise 

(dn(x)Sn(y)) = 2L 6n 2 S(x - y) (28) 

where the fluctuations are fully correlated for length scales less than Lq and Sn 2 is the average squared fluctu- 
ation. Simulations using these analytical solutions show that the system becomes depolarized for large fluctua- 
tions, the flip probability Pturb ~ 0.5. 

A more physically reasonable model of turblence is provided by the Kolmogorov spectrum in which the large 
scale density fluctuations cascade into density fluctuations on a small scale. These density fluctuations exhibit 
a power-law spectrum 
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dx(Sn(0)5n(x))e- lkx = C k a (29) 

where n(x) is the density of the supernova at the position x, k — 2it/\ where A is the wavelength of the 
fluctuation, Cq is a constant and a — —5/3. The density fluctuations on a large scale (the size of the shock 
wave) can be taken from numerical simulations and used to determine the value of Cq. Then this can be used 
to estimate the size of the fluctuations at the scale at which dominate the oscillations. Numerical solutions for 
the case of Kolmogorov turbulence also show that the system becomes depolarized for large turbulence [28] . 
Currently no exact solutions exist for Kolmogorov turbulence but the flip probability can be calculated in the 
perturbative limit giving 

! 0.84G F CWAm 2 sin2e w2/3 fnr ^(Am 2 sin2e) 2 
7T Am sin 20) 
1 f° r 1 IEW\ <<:1 

This is in the perturbative limit but a depolarization criterion can be defined by Pturb ~ 0.5. Fluctuations 
which satisfy this condition are expected to depolarize the neutrino flavours which are maximally mixed at the 
resonance. The current simulations show that this criteria is satisfied by a large margin for both the atmospheric 
and sterile resonances. In these cases the resonant neutrinos are expected to be depolarized in the turbulent 
region behind the shock wave. 

The flip probability of the density profile including the effect of turbulence is shown in Figs. [T]and[2j Note 
that for the forward shock the initial rise is unchanged by turbulence but the fall is reduced by the turbulent 
depolarisation. This just reflects the fact that the turbulence follows the shock wave. The effect on the reverse 
shock is dramatic, the depolarisation wiping out the reverse shock structure. Note too that the depolarisation 
following the shock front can lead to interesting observable effects in active neutrino oscillation for the case there 
are additional sterile neutrinos. This is because the active-sterile oscillation will destroy the near degeneracy 
between v & and v x . As a result the subsequent atmospheric neutrino oscillation signal, which vanishes if the 
active neutrino luminosities and spectra are equal, can be significantly enhanced. We will demonstrate these 
effects in the subsequent sections. 



3 Observing Supernova Neutrinos in Water 
3.1 Energy Spectra of Supernova Neutrinos 

In a type-II supernova approximately 3 x 10 53 ergs of energy is released and about 99% of this is in the form 
of neutrinos. The energy spectra of these neutrinos is determined by their interactions with matter. As the 
effect of weak magnetism in muon or tau production can be neglected, the spectra of v^v^v^ and v T are 
approximately equal and are collectively denoted by v x . The initial spectra of neutrino species v a from a 
supernova is parameterized as [42] 

, s S(t) (a(t) + lW'H 1 / E \ a{t) / E \ 

where (E) and $ are the average energy and total number flux and a is a dimcnsionlcss parameter which 
typically takes the values 2.5-5. For the results presented in this paper, we have assumed a g = 3 and a x = 4. 

To date only the Lawrence Livermore (LL) group [21] has published detailed results of the energy spectra of 
neutrinos for the duration of the supernova. In the LL simulations the luminosities of all flavours of neutrinos are 
approximately equal for times post bounce > 0.1 s. However there is a distinct difference in the average energies. 
This is due to the additional charged current interactions of the v e and v e with ambient matter. However these 
simulations did not include all the neutrino interactions now believed to be important |20j . Ongoing work 
by the Garching group including all possibly significant interactions predict the neutrino fluxes and average 
energies labeled Gl and G2 in Tabled This shows that the average energies between neutrino flavours become 
very similar and the luminosities are no longer equal. In order to compare the impact of the uncertainties on 
the average energies and fluxes of the neutrinos obtained in different supernova computer simulations, we will 
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Table 2: The average energies and total fluxes characterizing the primary neutrino spectra produced inside 
the supernova. The numbers obtained in the Lawrence Livermore simulations are denoted as LL, while those 
obtained by the Garching group are denoted as Gl and G2. 



present our results using supernova neutrino parameters given by both the Lawrence Livermore and Garching 
groups. Specifically, we consider the three cases shown in Table [2] [25]. For further details, we refer the reader 
to our earlier paper [15] . 

3.2 Signal in Water Cerenkov Detectors 

Supernova neutrinos with energy in the MeV regime will be dominantly detected in water through the capture 
of v e on protonsQ 

i> e + p -> n + e+ . (32) 

The emitted positron will be observed in the detector and its energy and time measured. Hence we should be 
able to get a fairly good reconstruction of the incoming v e energy spectrum and time profile. Neglecting the 
recoil of the neutron the energy of the neutrino E is related to the energy of the positron through the relation 
E e = E — 1.29, where all energies are in MeV. The number of events expected in a water Cerenkov detector 
from a galactic supernova explosion is given by 

N ,00 ,00 

N= — ±- / / F(E) * a(E) * e(E) * R(E - 1.29, E e ) dE dE e (33) 
4irD 2 J J 

where D is the distance of the supernova from earth, Nt is the number of target nucleons in 1 megaton of water, 
E is the energy of the neutrino, E e is the measured energy of the positron, F(E) is the flux at the detector, as 
defined in Eqs. (30) and (37), a(E) is the cross section, e(E) is the efficiency of detection and R(E — 1.29, E e ) 
is the energy resolution function. The efficiency is assumed to be perfect above 7 MeV and vanishing below this 
energy. The energy resolution function for which we assume a Gaussian form 

1 f —{Et — Eg)" 2 

= eX P r— J 

1tt<je V 2ct e 

where Et and E e are respectively the true and measured energy of the positron and we take the HWHM 



<je(Et) = VEqEt, where Eq = 0.22MeV. This is the same efficiency and energy resolution as used in [25| [43] . 
Water Cerenkov detectors usually are expected to have time resolution which is of the order of nanosecond. In 
what follows, we will bin our data either in time bins of 100 ms (at later times) or 10 ms (at earlier times). 
Since we expect the time resolution of the detector to be at least 3-4 orders of magnitude better, we do not 
include any time resolution function in our calculation of the number of events. 

In addition to the total number of events, we also calculate the the average energy of the detected positrons 
through the expression 

, F v = Jo" Jo" E A * HE) * a(E) * e(E) * R(E - 1.29, E A ) dE dE A 

{ e) f °° / °° F(E) * <t(E) *£(£)* R{E - 1.29, E A ) dE dE A [ ' 



4 Water Cerenkov detectors have other detection channels whereby they can observe u e and v x (electron scattering and charged 
and neutral current interactions on le O). However, the cross-section for these processes are much smaller and hence they are not 
considered here. 
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We show both the number of events as well as the average energy as a function of time. We also show 
the statistical uncertainties expected in 1 megaton water Cerenkov detectors. The statistical error in the total 
number of events are estimated as 



<r N = VN , (36) 
while that in the average energy is calculated as 



N 



(37) 



where ct/e) is the error in the average energy, N is the number of events, (E) is the average energy and (E 2 ) is 
the average energy squared. 

The neutrino flux in the detector is 

Fp = Y, F « P ^ ( 38 ) 

a 

where P a p is the oscillation probability and is given in Eqs. pg )) -(|25 ]) and F° is the initial flux of v a given in 
Eqn. (EH)- 



3.3 Input Supernova and Oscillation Parameters 

In what follows, we will present results for the typical values for the fluxes and average energies given in Table 
[2] for the LL, Gl and G2 "models". For the neutrino oscillation parameters, we assume the best-fit values 
Am^ = 8 x 1CT 5 eV 2 , sin 2 9i 2 — 0.31 and Am 2 ^ = 2.5 x 10~ 3 [ID] . For 9i 3 , we assume that it is large enough 
so that away from the shock, the Am^ driven resonant transition is fully adiabatic. Typically, this would be 
satisfied for sin 2 #i 3 10~ 3 . 

The other oscillation parameters relevant for supernova neutrino oscillations in the 3+2 scenario are con- 
strained by the combined data from Bugey, CHOOZ, CCFR84, CDHS, KARMEN, NOMAD, and LSND [29][44]. 
If we restrict the mass squared differences to lie in the sub-eV regime then the best-fit comes at Am^ = 0.46 
eV 2 , Amli = 0.89 eV 2 , U ei = 0.09, U e5 = 0.125, = 0.226 and = 0.16 [3l]. These values of the 
oscillation parameters give fully adiabatic transition at the resonance in the supernova when the shock is not 
present. 



3.3.1 Sterile neutrino sensitivity 

In the detailed estimates presented below we use the sterile neutrino parameters consistent with an explanation 
of the LSND experiment. However it is important to stress that the sensitivity to sterile neutrinos is much 
better than is needed to probe the LSND range and that, even if MiniBoone should rule out the sterile mixing 
angle regime used in our estimates, supernovae neutrino signals will still be important in the search for evidence 
for sterile neutrinos. It is straightforward to quantify the range of sensitivity. To observe the time dependent 
effects in the signal due to sterile neutrinos, the adiabaticity of the sterile resonances needs to be changed by 
the shock wave and/or the turbulence. For mixing angles sin 2 9ij <4x 10~ 6 the sterile resonances would be 
non- adiabatic for the entire time of interest. As a result there would be no oscillations into sterile neutrinos 
and the signal would be equivalent to that of only 3 active neutrinos. For sin 2 % > 4 x 10~ 4 the resonance 
would be adiabatic in the regions behind (with no turbulence) and in front of the shock wave. In this range the 
active-sterile resonant mixig effects are measureable. The situation is summarised in Table [3] 

For 4 x 10" 6 < sin 2 6 l3 < 4 x 10~ 4 there is a "transition region" , in which the effects discussed in this paper 
could be observed but may be less prominent. To quantify this first note that within this region of parameter 
space the mixing angles are sufficiently large such that each resonance is adiabatic in the absence of the shock 
wave. As well as changing the adiabaticity of each resonance the mixing angles change the relative proportion 
of each flavour eigenstate in each mass eigenstate, as a result changing the v e flux, F e in the detector. The 
approximate F e for small sterile mixing angles is shown in Table 4. If the sterile mixing angles and #13 are 
scaled as #y — > k&ij, where i=l,2, and j=3-5, and k < 1, F e is approximately unchanged except for the N2+I3 
and H2+N3 mass hierarchies. If the flux is unchanged both the number of events and the average energies are 
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sin 2 8ij 


Region 


6 x 10~ 3 < sin^ % 
4 x 1CT 4 < sin 2 By 
4 x 1(T 6 < sin 2 % < 4 x 1CT 4 
sin 2 0ij < 4 x 1CT 6 


Sensitivity range of MiniBoone 
Maximal effect of shock 
Transition region 
No effect of shock 



Table 3: The effect of the shock wave for the parameter space of sin 2 6^ with i=l,2 and j=3-5. 



unchanged. For the cases of N2+I3 and H2+N3 mass hierarchies F nos h oc k is scaled as F nos h oc k — ► k 2 F nos hock- As 
a result the flux and therefore number of observed events before and after the propagation of a shock wave scale 
as k 2 , the statistical uncertainties scale as A; -1 and the average energies remain unchanged. Therefore for k < 1 
the number of events decreases, and the uncertainties in the number of events and average energies increases. 
As the shock wave passes through the resonance the flux is F e ~ F s h oc k, which is independent of k. Therefore 
the total number of events is approximately unchanged. However, as the shock passes through a resonance the 
number of events increases, these corresponds to the lowest energy neutrinos, as a result the average energy 
decreases. For smaller k the relative increase in the number of events is larger and therefore the decrease in the 
average energy is larger. At later times there is an increase in the average energy as the resonance condition is 
satisfied for higher energy neutrinos, for smaller k the relative increase is larger and therefore the increase in the 
average energy is larger. As a result the structure of the average energy plot remains but is stretched for smaller 
k. Simulations show that the number of events during the shock propagation is approximately independent of 
k for k < 1 and the average energy plot is stretched as described above. 

3.4 Three active neutrinos 

We first consider the case of a "standard" supernova at 10 kpc from earth. If the mass hierarchy of the neutrinos 
is inverted, the flux of anti-electron neutrinos in the detector is given by 

F e = P° + P ee (P° - P°) (39) 

where P ee is given by Eq. (fl~5|) . For the values of sin 2 #i 3 ^, 10~ 3 that we assume throughout this paper, Pi 3 ~ 
in the absence of shock and the flux of v e at the detector F e — F®. As the shock crosses the resonance density, 
Pi 3 ~ 1 as discussed in section 2.2 and the flux of D e at the detector is given by F e ~ (1 — \U e \\ 2 )F x + \U e i\ 2 F° . 
If the average energy of v x is larger than that of v e as obtained in the LL simulations, the average energy will 
decrease as a result of the shock effect. This can be seen from Fig. [3^a), where we show the time evolution of the 
average energy of the neutrinos detected in a megaton water Cerenkov detector. The band shows the statistical 
error expected in the measured average energy. For no shock (NS) the average energy remains almost constant 
with time. On the other hand the average energy decreases and hence shows "bumps" in the time profile as the 
shock crosses the position of the Amf^ driven resonance. We get a single bump for the forward shock only and 
double bump when the reverse shock is also present. This is due to the shape of the flip probability shown in 
Figs. [T]and[2] Also shown (in Fig. [3(c) and Fig. [3je)) is what is expected if the difference between the average 
energies of the T> e and v x created in the supernova were not so well separated. We notice that even though the 
effect of the shock is less dramatic, nonetheless it is still there and should be observable. Particularly note that 
for G2 the average energy of v e and v x are equal, therefore the change in the average energy measured by the 
detector is due to the change in the number of detected neutrinos. In Fig. 2]we show the total number of events 
expected for the LL (panel (a)), Gl (panel (c)) and G2 (panel (e)) cases, with and without the shock effects. 
Since the cross section for the detection cross section increases quadratically with energy, the number of events 
will decrease when the shock crosses the Am§i resonance density and this results in lowering of the number of 
events and hence again gives single bump and double bump for the forward and forward+reverse shock cases 
respectively, as a function of time. 

The corresponding results taking into account the effects of turbulence are shown in right-hand panels of 
Figs. [3jand[4] Panels (b), (d) and (f) in these figures correspond to the LL, Gl and G2 simulations respectively. 
The inclusion of turbulence changes the signal for t > 5s. This is when the resonant density is in the turbulent 
region. As discussed before, the system is largely depolarized giving Pi 3 ~ 1/2. The flux at the detector when 
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Table 4: The flux of neutrinos in the approximation that |C/ e i| 2 , | C/ e 2 1 2 >> | t/ e 3 1 2 , \U e 4,\ 2 , \U e $\ 2 and Pu,P3i — 
P 2 i, where i=4 or 5. The total flux in presence of shock is given by F e = F noshock + F shock , where F noshock = F e 
without the shock and F s hock is the extra component due to the shock effect. 

the shock crosses the resonance region is now given by F e ~ (1 — 0.5\U e i\ 2 )F x ) + 0.5|C/ e i | 2 F e °. For the same 
arguments as before, the average energy is lower than the case of no shock, however the effect now is much less 
than the case where turbulence was not considered. Therefore, the single and double bump features expected 
from shock effects get smeared out to a large extent. However, we still see a non-trivial variation of the expected 
average energy in the detector as a function of time, both for the LL and G1/G2 simulations. 

3.5 Three active and two sterile neutrinos 

For the case where we have two extra light sterile neutrinos in addition to the three active ones, the situation 
gets a lot more involved due to the possibility of multiple resonances as discussed in section 2.1. The effect of 
the shock is also richer here since the shock passes through the different resonance densities at different times. 
In Table 3] the resultant neutrino flux at the detector is given in terms of the original neutrino fluxes produced, 
for the different neutrino mass spectra considered, for both with and without shock effects. We reiterate that 
we have chosen mixing angles such that in absence of shock effects all resonances are adiabatic. The effect of 
the shock is to turn the adiabatic resonant transition into non-adiabatic ones, as discussed before. For the case 
where sterile neutrinos are present, the change in number of events are characterized by the oscillations of active 
neutrinos into sterile species and vice-versa. Thus in this case, difference in the initial neutrino energy spectra is 
not a pre-requisite for observing resonant oscillations and shock effects, unlike in the case of 3 active neutrinos 
only. As the shock moves in time, its effect is imprinted on the time dependence of the neutrino signal, as in 
the case of 3 active neutrinos only. However, with sterile neutrinos there are further modulations in the signal 
because the shock passes through the additional sterile resonances. The number of resonances is dependent on 
the number of sterile neutrinos as well as the neutrino mass spectrum and therefore the expected signal is often 
different for each scenario, as described in detail in [19 . In fact, as we will see, even when the shock passes 
through the Am^ driven active resonance, the time profile and energy spectrum of the observed neutrinos in 
many of the possible mass spectrum listed in Eqs. (|17p - (j2"2")l , are different than what is expected for the case of 
3 active neutrinos only. Typically, the shock crosses the multiple "sterile resonances" at very early times (t & 2 
sec), while it crosses the Am^ driven "active resonance" at later times (t ^ 3 sec). Therefore in what follows, 
we will present results separately at late and early times to show clearly the effect of the shock wave on the 
neutrino signal in the detector. For later times, we will consider time bins of 100 ms, while at earlier times since 
the time dependence is much more sharp, we present our results for smaller time bins of 10 ms. 

3.6 Neutrino signal at late times 

We begin by discussing the evolution of the expected neutrino event rate and average energy at later times. The 
total number of events for the different mass spectra are shown in Fig. [5] while the expected average energies 
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are shown in Fig. at times between t = 3 — 12s. These plots show the impact on the signal when the shock 
wave passes through the Am^ resonance between active neutrinos. In panels (a)-(d) of Figs. [5] and [5] we show 
only the mass spectrum cases where Am^ < 0, in which we have resonance and hence also shock effects. Which 
mass spectra will get the shock effects can be easily seen from Table |4j Since at late times the resonance that 
gets affected by the shock is the Ato^ driven resonance between the 1-3 states, the relevant jump probability 
involved is P13. All other jump probabilities are zero here, from our choice of the mixing angles. We can see 
from TableUJthat the N2+I3 and 12+13 are the only cases which will get affected. For all the other mass spectra 
possible, we do not expect any modulation in the signal at late times due to shock effects. Hence we show only 
the signal for the N2+I3 and 12+13 cases. For comparison we show in the last 2 panels of these figures the case 
of only active neutrinos (13). 

In the approximation that |t/ e 3| 2 , |£/ e 4| 2 and If/^ 2 can be neglected in comparison to | £/ e i | 2 and |t/ e 2| 2 , we 
note from Table 2] that in absence of shock, 



13 F e ~ F° , (40) 

N2 + 13 => F e ~ , (41) 

12 + 13 => F e ~ F° , (42) 

whereas when the shock passes through the 1-3 resonance the fluxes are modified to, 

73 =► F e ~ F°(l - P 13 \U el \ 2 ) + P 13 |f/ el | 2 F e ° , (43) 

N2 + I3 =► F e ~|[/ e i| 2 P 13 F° , (44) 

12 + 13 F e ~ (I - P 13 \U el \ 2 )F° . (45) 



By comparing Eqs. (|40[) - (|42p we note that in absence of shock, while the fluxes are same for 13 and 12+13 
spectra, N2+I3 predicts almost zero fluxes. The black bands in Fig. [5] corroborate the above statement. When 
the shock passes through the Am^ resonance it increases the P e flux for N2+I3, while it decreases the same for 
12+13, as can be seen by comparing Eqs. (|4"Tj) and (|4*2"|) with Eqs. (|4"4"|) and (|4"5|) . In the case of 13, the resultant 
flux is an admixture of a fraction of the initial v e and v x fluxes. However, since the average energy of initial v e 
flux was smaller, the total number of events goes down in the detector when the shock passes through the 1-3 
resonance even in this case. We stress that the sterile cases are qualitatively different from the 13 case, since for 
them the net number flux sees a big increase or decrease due to shock. Therefore, while for 13, the effect of the 
shock wave comes predominantly through the difference in the average energy of v e and u x , for sterile cases we 
see a combined effect coming from a direct change in the number of neutrinos arriving on earth as well as the 
difference in the energy spectra of the different species. 

For all the mass spectra shown, under the assumption that there is no turbulence, typically a single bump 
is observed for a forward shock and a double bump for a forward and reverse shock, as expected from the shape 
of the flip probability shown in Figs. [T]and[2J The effect of taking the turbulence into account is to smear these 
sharp changes in the oscillation probability due to the effective depolarizing of the resonances. With turbulence, 
a single bump is observed, followed by a region in which the number of events are typically different from what 
is expected for no turbulence. This can be seen in the right-hand panels in Fig. \5\ 

The expected average energies for the N2+I3 and 12+13 cases also show a very striking evolution with time 
at < ^ 3s, which is very different from that predicted by 13. In the case of 13, the average energy decreases in 
presence of shock as expected, since the shock reduces the conversion probability of the high energy v x into D e , 
thereby decreasing the average energy. Thus it usually shows 1 sharp decrease for forward shock and 2 sharp 
falls for forward+reverse shock. On the other hand both N2+I3 and 12+13 predict that the average energy 
fluctuates on both sides of the average energy expected in absence of shock. The key issue to note here is that 
the position of resonance is determined by the energy of the neutrino. The higher (lower) energy neutrinos go 
through the resonance at lower (higher) density. Since the shock moves from higher to lower densities in time, 
the lower energy neutrinos are affected by the shock earlier than the higher energy ones. For the sterile cases, 
the effect is more subtle. Here the time evolution of the average energy is a combined effect of the change in 
the total flux as well as the energy dependence of the resonance position. For the 12+13 case since resonance 
happens for lowest energy neutrinos first, the effect of the shock is to reduce them in the flux (cf. Eq. (|42|) and 
(|45|) ). thereby increasing the average energy. Eventually, the shock goes through higher energy resonances, and 
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this then reduces the average energy. For the N2+I3 case also the resonance happens for the lowest energy 
neutrinos first, but now the shock effect increases them in the flux and thereby decreases the average energy. 
Eventually, the shock goes through the high energy resonances, increasing the average energy. 

The right-hand panels of Fig. [6] show the average energy after turbulence is taken into consideration. The 
presence of turbulence has a typical effect on the time evolution of both the total number of events and the 
average energy deposited by the (anti)neutrinos. One can see that without turbulence both the average energy 
and number of events after the passage of (both) shock(s) go back to its pre-shock value. This is because the 
P ee before and after the shock are exactly the same in this case. However, since the turbulence changes the 
flip probability permanently to Py behind the shock, the number of events and average energy are consistently 
lower once the shock crosses the resonance point. 

3.7 Neutrino signal at early times 

In Figs. [7] and [9] we show the number of predicted events and average energies as a function of time (in absence 
of turbulence) for the 5 different mass spectra, between t = 0.1 — 1.5s, in short time bins of 10 ms. Figs. [8] and 
[Tm show the corresponding results when turbulence is taken into account. We do not show results for the 12+13 
spectra here since, as can be seen from Table 21 it depends only on Pi 3 which will have shock effects only at 
the Am§! resonance at late times. 

As in the previous subsection, Fig. [Jjcan be understood in terms of the flux predictions in the presence and 
absence of the shock, given in Table [U The only difference is that at early times the shock passes through the 
sterile resonances and hence we expect contributions coming from the jump probabilities associated with the 
sterile resonances. In particular, we note that in the absence of shock, the predicted flux on earth is almost 
zero for N2+I3 and H2+N3, while for H2+I3 and I2+N3 we expect the flux to be \U e2 \ 2 F° and for N2+N3 it 
is predicted to be |^7 e i | 2 P^- Since \U e i\ 2 > \U e 2\ 2 , the expected signal is larger for N2+N3. These features are 
evident from Fig. UJ As a result of the shock we see a modulation in the resultant signal visible as bumps in the 
figure. Note that in all the 5 cases shown, the shock effect increases the number of events. The shock powered 
modulations get affected when turbulence is taken into account and the corresponding results are shown in Fig. 

m 

Fig. [9] shows how the average energy evolves at early times. As we had seen at later times, the effect of the 
shock wave is to make the average energy fluctuate on both sides of the corresponding static density case. The 
reason for the initial average energy decrease and subsequent increase is also the same as that discussed in the 
previous subsection. Fig. rrU]shows the average energy evolution when the effect of turbulence is considered. The 
rapid fluctuations are mellowed due to the presence of turbulence, however we still see statistically significant 
fluctuations in the average energy at very early times, a feature which, if observed, would provide an almost 
model independent signal of the presence of extra sterile neutrino species which are mixed mildly with the active 
neutrinos. 

As the water Cerenkov detector can measure the energy of the incoming neutrino rather efficiently, the events 
can be binned in energy as well as time. In Fig. [IT] we present the number of events expected in energy bins 
of width 10 MeV and time bins of 10 ms. We show results for the N2+I3 spectrum only as an example. This 
figure is shown for only early times to illustrate the effect of the sterile resonances. As the resonance densities 
are energy dependent, the propagation of the shock wave can be observed when the events are binned in energy 
and time. The shock wave crosses the resonant density corresponding to neutrinos with lower energies first, 
therefore the characteristic 'bumps' are observed in the lower energy bins first. The presence of such energy 
and time dependent bumps in the resultant signal in the detector provides a 'smoking gun' signal for sterile 
neutrinos. 

4 Comparison and Discussions 

In Table [5] we show the model independent characteristic features in the expected signal for the eight different 
cases considered in this paper, the two only active cases (N3 and 13) and the six active plus sterile cases. The 
time interval is divided into early (t ^ 2 s) and late (i <; 3 s) times, as before. For each case we state if we 
expect large (L) or very small (S) number of events in the 2 time zones. We also indicate in the table whether 
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Table 5: The expected trends in the neutrino signal in a megaton water Cerenkov detector at early and late 
times. The symbols L and S signify large and small number of events respectively. The presence or absence of 
shock effects are shown by the / and x symbols. 

we expect shock induced sharp modulations in the signal or not for each of the possible mass spectra. By 
comparing the predictions for all the cases, we can infer the following: 

• 1. If large events are seen at early times, 

2. shock effects arc not seen at early times, 

3. a large number of events are seen at late times and 

4. shock effects are not seen at late times 
then we must have the N3 spectrum. 

• 1. If very few events are seen at early times, 

2. shock effects are seen at early times, 

3. very few events are seen at late times and 

4. shock effects are seen at late times 
then we must have the N2+I3 spectrum. 

• 1. If very few events are seen at early times, 

2. shock effects are seen at early times, 

3. very few events are seen at late times and 

4. shock effects arc not seen at late times 
then we must have the H2+N3 spectrum. 

• 1. If a large number of events are seen at early times, 

2. shock effects are seen at early times, 

3. a large number of events are seen at late times and 

4. shock effects are not seen at late times 

then we could have either N2+N3 or H2+I3 or I2+N3 spectrum. 

• 1. If a large number of events are seen at early times, 

2. shock effects are not seen at early times, 

3. a large number of events are seen at late times and 
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4. shock effects are seen at late times 

then we could have either 13 or 12+13 spectrum. This degeneracy can be split by measurng the energy 
dependence (c.f. Figure 6 (c) and (d)). 

Therefore, from the Table [5] we conclude that, irrespective of model uncertainties, the presence of sterile 
neutrinos can be easily proved from modulation of the signal due to shock effects at early time. Only the 12+13 
case does not give any time modulation, even though the sterile neutrinos are present in this case. In addition, 
from the point-wise discussion above we can see how well we can distinguish one mass spectrum from the other. 
N3, N2+I3 and H2+N3 can be uniquely determined by comparing early time behavior of the signal with its 
late time behavior. The N2+N3, H2+I3 and I2+N3 can be separated from the rest, but since they predict 
similar trends at early and late times, one will be need a more careful model dependent study to unambiguously 
disentangle them from each other. Similarly, though 13 and 12+13 can be separated from the other spectra, one 
will need a more careful analysis to distinguish between the two. 

5 Summary and Conclusions 

Model independent information about the neutrino mass spectrum can be obtained through observation of 
signatures of supernova shock wave(s) on the resultant neutrino signal in terrestrial detectors. In particular, 
such signals probe the existence of extra sterile neutrino through their possible resonant transition with active 
neutrinos inside the supernova. Oscillations between sterile neutrino and active ones are characterized by unique 
signatures in the final neutrino energy spectrum as well as their evolution with time. While the impact of the 
presence of sterile neutrinos on the time evolution of the neutrino is due to shock effects, their impact on the 
resultant neutrino flux and spectra comes from neutrino oscillations both with and without shock effects. 

Concerns had been raised recently over the observability of the shock effects due to the turbulent density 
variations following the shock wave which prove to be very important. In this paper we made a detailed study 
of the turbulent shock effects on the neutrino induced galactic supernova signal expected in megaton water 
Cerenkov detector. Water detectors can give information on the number, the energy spectrum, as well the time 
profile of the arriving (anti)neutrinos. We have shown that the impact of the supernova shock waves is evident 
in the neutrino signature in megaton water Cerenkov detectors in all of these, making this class of detectors 
extremely good for studying shock effects. We have considered the impact of the turbulence left behind by the 
shock wave and have seen that although the shock effect is diluted, it is still significant and we still expect to 
observe them in the megaton class of detectors. We have concentrated on the discernable signatures of sterile 
neutrinos, which might be mixed mildly with the 3 active neutrinos. We have illustrated the effects using sterile 
neutrino parameters which were chosen in a fit to all neutrino data including LSND. However we showed that the 
observable signals persist for much smaller mixing angles than are observable by the LSND (or MiniBOONE) 
experiments. Hence our results are relevant whatever the outcome of the MiniBOONE data. 

The most striking evidence for sterile neutrinos in the supernova neutrino signal are sharp bumps at t ^ lsecin 
the observed number flux as well as the average energy of the v e detected through their capture on protons in 
megaton water detectors. These can be caused only by sterile neutrino resonances inside the supernova. Only 
the 12+13 case for the mass spectrum does not predict these early time shock induced modulations in the signal. 
In addition, a model independent comparison of the signal trend between early and late times can give us a 
rather unambiguous signature on 3 of the 8 possible mass spectra considered in this paper. The N3, N2+I3 
and H2+N3 cases predict unique combination of behavior at early and late times and hence can be determined 
model independently from the observations. The remaining 5 cases can be classed into 2 categories depending 
on their combination of predicted trends at early and late times. Distinguishing the 13 from the 12+13 spectra 
would require a more careful model dependent analysis of the future supernova neutrino data. Similarly, the 
N2+N3, H2+I3 and I2+N3 can be separated from the relative differences in their predictions of average energy 
and number of events as a function of time. 

One important byproduct of the turbulent effects involving active-sterile neutrino mixing is that, even if the 
energy spectrum and luminosisites of the active neutrinos are initially the same, the depolarising effect of the 
turbulence for the active sterile resonances in the wake of the shock front will make the active neutrino spectra 
spectra and luminosities significantly different. As a result the atmospheric neutrino resonant effects involving 
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the active neutrinos in the presence of sterile neutrinos may be expected to give rise to more significant effects 
in the supernova neutrino signal than is the case without sterile neutrinos. 
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Figure 3: The average energy in 100ms bins, for the cases of a forward shock (FS), a forward and reverse shock 
(FRS) and no shock (NS): (a)LL no turbulence, (b)LL with Turbulence, (c)Gl no turbulence, (d) Gl with 
turbulence, (e)G2 no turbulence, (f) G2 with turbulence. 
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Figure 4: The number of events in 100ms bins, for the cases of a forward shock (FS), a forward and reverse 
shock (FRS) and no shock (NS): (a)LL no turbulence, (b)LL with Turbulence, (c)Gl no turbulence, (d) Gl 
with turbulence, (e)G2 no turbulence, (f) G2 with turbulence. 



21 



(a) 



^ Gl N2+I3 FS no Turbulence 
I Gl N2+I3 FRS No Turbulence 
fl^B Gl N2+I3 MS no Turbulence 




2500 



c 2000 



^ Gl 112+13 FS Turbulence 
| 51 K2+I3 FRS Turbulence 
I 51 K2+E3 us Turbulence 




Time (s) 



rime (s) 



□ 


;;i 


12+13 


FS Ho Turbulence 




si 


12+13 


FRS Ko Turbulence 




;;i 


12+13 


ns no Turbulence 




2000 



1000 



a 5 ny 



^| Gl 12+13 FS Turbulence 

Gl 12+13 FRS Turbulence 
■ 51 I2+I3 hS Turbulence 




Time (s) 



rime (s) 



(e) 



I I Gl 13 FS mo Turbulence 
| Gl 13 FRS No Turbulence 
■ Gl 13 ns no Turbulence 




10 



2500 



I I Gl 1 3 FS Turbulence 
| Gl 13 FRS Turbulence 
H Gl 13 ns Turbulence 



Time (s) 



(a a 
Time (s) 



10 




Figure 5: The number of events in 100ms bins, for the cases of a forward shock (FS), a forward and reverse shock 
(FRS) and no shock (NS): (a)N2 + 13 no turbulence, (b)N2 + 13 with turbulence, (c)I2 + 13 no turbulence, 
(d)I2 + 13 with turbulence, (e)I3 no turbulence, (f)I3 with turbulence. 
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Figure 6: The average energy in 100 ms bins, for the cases of a forward shock (FS), a forward and reverse shock 
(FRS) and no shock (NS), with the initial energy spectra Gl as discussed in the text: (a)N2 + 13 no turbulence, 
(b)N2 + 13 with turbulence, (c)I2 + 13 no turbulence, (d)I2 + 13 with turbulence, (e)I3 no turbulence, (f)I3 
with turbulence. 
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Figure 7: The number of events in 10 ms bins, for the cases of Amfj = 0A6MeV 2 and Am^ = 0.54MeV 2 , 
with no turbulence, a forward shock (FS), a forward and reverse shock (FRS) and no shock (NS), with the 
initial energy spectra Gl as discussed in the text:(a)N2 + N3, (b)N2 + 13, (c)H2 + N3, (d)H2 + 13, (c)I2 + 
N3. 
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Figure 8: The number of events in 10 ms bins, for the cases of Am^ = 0A6MeV 2 and Am| x = 0.54MeV 2 , 
with turbulence, a forward shock (FS), a forward and reverse shock (FRS) and no shock (NS), with the initial 
energy spectra Gl as discussed in the text:(a)N2 + N3, (b)N2 + 13, (c)H2 + N3, (d)H2 + 13, (c)I2 + N3. 
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Figure 9: The average energy in 10 ms bins, for the cases of Am^ = 0A6MeV 2 and Am^ = 0.54MeV 2 , with 
no turbulence, a forward shock (FS), a forward and reverse shock (FRS) and no shock (NS), with the initial 
energy spectra Gl as discussed in the text:(a)N2 + N3, (b)N2 + 13, (c)H2 + N3, (d)H2 + 13, (e)I2 + N3. 
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Figure 10: The average energy in 10 ms bins, for the cases of Am^ = 0A6MeV 2 and Arn^j = 0.54MeV 2 , with 
turbulence, a forward shock (FS), a forward and reverse shock (FRS) and no shock (NS), with the initial energy 
spectra Gl as discussed in the text:(a)N2 + N3, (b)N2 + 13, (c)H2 + N3, (d)H2 + 13, (e)I2 + N3. 
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Figure 11: The number of events in 10 ms and 10 MeV bins, for N2 +13 and the cases of a forward shock 
(FS), a forward and reverse shock (FRS), with the initial energy spectra LL as discussed in the text:(a)FS no 
turbulence, (b)FRS no turbulence, (c)FS with turbulence, (d)FRS with turbulence. 
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